Molecular Oncology

BCL2L13 is a mitochondrial BCL2 family protein linked to mitophagy and ceramide metabolism, but its role in NSCLC metastatic plasticity remains unclear. Human lung cancer Tissue Microarray and matched patient specimens showed subtype and site dependent BCL2L13 expression, with higher cytoplasmic granular staining in primary NSCLC and reduced, heterogeneous staining in lymph node metastases, most evident in adenocarcinoma and squamous cell carcinoma. Because Epithelial mesenchymal transition and anoikis resistance are central requirements for metastatic dissemination, this primary to node attenuation provided the rationale to test BCL2L13 knockdown and overexpression in metastasis relevant NSCLC models. In A549 and LLC cell lines. TGF beta 1 induced coordinated mitophagy and EMT with mitochondrial enrichment of BCL2L13. BCL2L13 knockdown impaired TGF beta 1 and carbonyl cyanide m chlorophenyl hydrazone associated mitophagy, reducing LC3 beta mitochondria colocalization, TOMM20, LAMP1 overlap and mitochondrial LC3 II, p62, TOMM20 turnover; BNIP3 and NIX redistribution did not compensate. BCL2L13 loss enhanced EMT marker switching and migration, whereas overexpression partially opposed these changes. During detachment, BCL2L13 knockdown reduced anoikis associated apoptosis despite preserved mitochondrial recruitment of BAX, BAK, BNIP3,NIX, altered BID processing, non parallel caspase activity and shifted FAK phosphorylation. Pharmacological autophagy modulation did not reverse this anoikis phenotype. Lipidomics identified adhesion state dependent ceramide synthases CerS2, CerS6 linked sphingolipid remodeling: BCL2L13 knockdown increased C24 linked sphingolipid species in attached cells but reduced C16, C24 ceramide related profiles during anoikis. These findings identify BCL2L13 downregulation as a metastasis associated mitochondrial-lipid state that limits mitophagic quality control while favoring EMT and detachment survival in NSCLC adenocarcinoma.

Background: Spatial organization is increasingly recognized as a key determinant of tumor-immune interactions in head and neck squamous cell carcinoma (HNSCC). The GSE300147 Xenium spatial transcriptomic resource generated by McCord and colleagues established a framework for mapping spatially coordinated T-cell states in HNSCC. However, how tumor-enriched epithelial immune states relate to metabolic, redox, and stress-adaptive transcript programs remains incompletely defined. Methods: A secondary, data-driven reanalysis of GSE300147 was performed, focusing on 17 confirmed HNSCC Xenium sections after exclusion of a non-HNSCC ameloblastoma specimen. A total of 1,148,244 cells were analyzed, including 558,867 EpCAM+ tumor-enriched epithelial cells. Tumor-enriched epithelial cells were classified into Hot, Intermediate, and Cold states using a Composite Hotness framework integrating T-cell inflammatory signature score, checkpoint-associated signaling, CD274 expression, IFN/antigen-presentation signature score (IFN/AP), and tumor-immune proximity. Six metabolic ecosystem states, neighborhood profiling, spatial permutation testing, and an integrated Immune-Metabolic-Redox Ecosystem Score (IMRES) were then applied. Results: Immune activation was spatially heterogeneous across HNSCC sections. Immune-hot tumor-enriched epithelial regions showed not only inflammatory, checkpoint-associated, and antigen-presentation signature scores, but also coordinated metabolic, oxidative-redox, and stress-response transcript programs. IMRES, derived from available immune, metabolic, redox, and stress-response transcript components represented in the Xenium panel, increased progressively from Cold to Intermediate to Hot tumor-enriched epithelial states and was associated with NFE2L2, GDF15, HLA-DRA, CD274, KEAP1, and MDM2. Integrating IMRES with Composite Hotness identified a distinct Hot+IMREShigh ecosystem comprising 106,874 tumor-enriched epithelial cells. This state showed the strongest immune-active and stress-adaptive features and was positioned closer to immune populations than expected by random assignment. An alternative rank-based robustness analysis reproduced the IMRES-associated ecosystem axis and correlated with the original module-based score (Spearman r = 0.597). Conclusions: This secondary reanalysis extends the original spatial T-cell framework by defining a complementary tumor-centered immune-metabolic-redox ecosystem in HNSCC. IMRES provides a transcript-derived framework for identifying Hot+IMREShigh neighborhoods where immune activation, checkpoint signaling, metabolic remodeling, and stress adaptation converge, providing a hypothesis-generating framework for studying immune resistance and therapeutic vulnerability.

Bulk telomere length measured from tumour sequencing is routinely interpreted as a property of the cancer cells. However, a tumour specimen is a mixture, and the patient who supplies it has a telomere length of their own. Here I re-analyse published pan-cancer telomere estimates and ask how much of a tumour's telomere length is patient-specific. A calibration step comes first. Whole-genome and low-pass estimates recover the known cross-sectional attrition of leukocyte telomeres with age, at 26.6 bp per year in blood normals, whereas whole-exome estimates do not. After adjustment for cancer type, sequencing centre and sex, the exome slope is minus 0.6 bp per year. In 684 blood-normal aliquots sequenced by both assays, the whole-genome estimate declines at 38.9 bp per year, whereas the exome estimate from the same DNA shows no detectable decline. The difference between assays is 41.5 bp per year, with P = 3 x 10^-10. Because exome data constitute 78.6% of the original resource, downstream analyses use only whole-genome and low-pass libraries. Within those data, tumour telomere length tracks the patient's matched-normal telomere length. The Spearman correlation is 0.395 in TCGA, with positive associations in 22 of 23 cancer types. This finding replicates in PCAWG using a different telomere estimator, with a correlation of 0.472 and positive associations in all 24 histologies examined. Adjustment for cancer type, sequencing centre and library type leaves a regression coefficient of 0.385. The association is also stable after adjustment for age, sex, tumour purity, leukocyte fraction, ploidy, sequencing coverage and continental ancestry, with coefficients ranging from 0.406 to 0.429. Pure normal-cell admixture is rejected as the sole explanation. Under a two-compartment mixture model, the coefficient for host telomere length is expected to equal 1 and the host-by-purity interaction to equal minus 1. These restrictions are jointly rejected with P = 0.001. Tumour purity, leukocyte fraction and age each explain only about 1 to 3% of within-cohort variance and do not alter the cross-cancer ranking. By contrast, the between-cohort coefficient is not directly interpretable. Its apparent near one-to-one relationship with tissue-associated telomere length depends strongly on which tissue supplies the matched-normal reference and on the statistical spread of that predictor, falling to 0.44 when organ-matched solid tissue is used. Bulk tumour telomere length is therefore a composite phenotype containing a replicated patient-specific component. Telomere biomarker studies should include matched-normal telomere length as a covariate rather than treating tumour telomere length as exclusively tumour-intrinsic.

Large genomically annotated acute myeloid leukaemia (AML) datasets exist, but population-based contemporary cohorts remain scarce. Here we report clinicopathological, genomic, and outcome data from Danish AML patients. 2,512 AML patients were identified between 2015-2022, of whom 33.8% had available NGS data (NGS+). In patients [≤]70 years, baseline characteristics and outcomes were comparable between NGS+ and NGS- groups. In patients >70 years, more NGS+ patients received intensive treatment, but survival was similar among intensively treated patients. The distribution of mutations varied significantly by age and sex, with older age and male sex exhibiting higher frequencies of adverse-risk gene mutations. In intensively treated NGS+ patients, ELN2017 stratified 5-year OS: 58.4% (favorable), 43.4% (intermediate), and 28.2% (adverse), with hazard ratios (HRs) of 0.63 (favorable) and 1.45 (adverse) relative to intermediate. ELN2022 yielded corresponding OS rates of 56.9%, 51.8%, and 29.7%, with HRs of 0.78 and 1.86. The two models had comparable predictive performance for OS in a time-dependent model. In conclusion, outcomes of intensively treated AML patients were comparable irrespective of NGS status, underscoring the representativeness of the REFORM-AML database for the Danish AML population. Age and male sex correlated with adverse-risk mutations, and both ELN2017 and ELN2022 robustly predicted survival.

Abstract Background: Cytolytic activity (CYT), a widely used transcriptomic surrogate of anti-tumor immune cytotoxicity derived from GZMA (granzyme A) and PRF1 (Perforin 1) expression, is associated with clinical outcomes across cancers. MicroRNAs (miRNAs) are key post-transcriptional regulators of tumor immunity, yet their pan-cancer roles in modulating cytolytic activity remain incompletely understood. Objective: This study aimed to identify conserved miRNA regulators of tumor cytolytic activity and their downstream gene-mediated networks across diverse cancer types, while evaluating their clinical and therapeutic relevance. Methods: Matched miRNA and mRNA expression profiles from 9,288 primary tumors across 31 TCGA cancer types were analyzed. A multi-stage framework was applied: per-cancer Spearman correlations (|{rho}| >= 0.30, FDR < 0.05) identified recurrent CYT-associated miRNAs (at least 3 cancer types); these were integrated with TargetScan-predicted targets and subjected to pan-cancer and cross-cancer triple filtering (miRNA-gene and gene-CYT associations). All associations underwent tumor purity adjustment using Consensus Purity Estimate (CPE), with LUMP (Leukocytes Unmethylation for Purity) as sensitivity analysis. Candidates were further prioritized by random forest modeling with bootstrap stability, cancer-type-adjusted Cox regression, mediation analysis, immune cell deconvolution, k-means molecular subtyping, pathway enrichment, and DGIdb-based drug-target prioritization. Results: The analysis converged on 38 high-confidence miRNA-gene-CYT regulatory triplets involving 9 conserved miRNAs and 31 target genes after stringent purity adjustment and multi-layer validation. All nine miRNAs exhibited complete bootstrap stability. Mediation analysis confirmed significant gene-level mediation in 37 of 38 triplets (FDR < 0.01), with mediated proportions up to 94%. The final miRNA signature defined two distinct pan-cancer immune subtypes (immune-hot vs. immune-cold) with significantly different cytolytic activity and overall survival (OS) (HR = 0.754, FDR = 1.12 x 10^-4). The network was enriched for T-cell activation and lymphocyte differentiation pathways and highlighted multiple druggable targets, including CTLA4 and CD274 (PD-L1), nominating 124 candidate compounds. Conclusions: In conclusion, this tumor purity-adjusted pan-cancer study defines a compact, reproducible, and clinically relevant miRNA network that regulates cytolytic activity across diverse malignancies. By linking miRNA biology to immune subtyping and actionable therapeutic targets, the present work provides a valuable foundation for advancing precision immuno-oncology.

Lung squamous cell carcinoma (LUSC) is an aggressive malignancy characterized by high cellular plasticity and few targeted treatment options. TNIK overexpression is common in LUSC and promotes tumor growth, with TNIK inhibition sensitizing LUSC to radiotherapy, though the underlying mechanisms are not well defined. Through transcriptomic analyses and functional assays, we identified TNIK as a regulator of a MYC-dependent transcriptional network that coordinates epithelial-mesenchymal plasticity and cell proliferation in LUSC. Depletion of TNIK reprogrammed LUSC cells from a hybrid epithelial/mesenchymal state towards an epithelial, senescent-like state characterized by reduced cell migration, invasion, reduced DNA synthesis, and enhanced {beta}-galactosidase activity. Using a small-molecule screen approach, we found that TNIK inhibitors cooperated with agents suppressing the histone methyltransferase and MYC binding partner EZH2, which further suppressed partial epithelial-to-mesenchymal transition (pEMT). Mechanistically, we identified MYC as a key downstream TNIK effector in LUSC cells: MYC depletion phenocopied the effects of TNIK loss on pEMT and senescence, and restoring MYC expression bypassed the effects of TNIK depletion. Collectively, these results implicate TNIK in the mechanisms linking epithelial-mesenchymal plasticity with proliferation and evasion of senescence and provide insights into future strategies for the clinical deployment of TNIK inhibitors in LUSC and other TNIK-dependent malignancies.

Introduction: Gliomas integrate into neural circuits and heighten neuronal excitability, engaging in bidirectional communication whereby neuronal activity promotes tumor growth and proliferation. Aging reshapes the brain microenvironment through extracellular matrix changes, altered secretory factors, and immune dysfunction, creating conditions permissive to tumorigenesis and limiting immunotherapy efficacy in glioblastoma. However, its effect on neuronal excitability and signaling in glioblastoma remains poorly understood. Methods: We developed a novel classification system for glioblastoma by leveraging three classes of DNA methylation-based aging biomarkers: chronological, biological, and mitotic clocks. This approach stratified tumors into accelerated and decelerated epigenetic aging subtypes, which we then characterized at the molecular, functional, and clinical levels using multimodal analyses. Guided by these profiles, we evaluated the in vitro effects of the FDA-approved agents levetiracetam and riluzole, alone and in combination with temozolomide, on U87MG and A172 cell lines. Specifically, we assessed changes in cell viability, apoptosis, and the expression of marker genes related to stemness, neuronal hyperexcitability, and immunosuppression. Results: Tumors with decelerated epigenetic aging showed expression modules and CpG hypomethylation associated with neuronal activity and stemness, and carried significantly worse prognosis. Single-cell and spatial multi-omics analyses revealed enrichment for neurons and malignant neural stem-like cells in these tumors. They also displayed enhanced intercellular communication, driven predominantly by glutamate signaling across the malignant, neuronal, and immune compartments of the tumor microenvironment. In vitro pharmacological inhibition of glutamatergic signaling with levetiracetam and riluzole reduced cell viability, induced apoptosis, and suppressed expression of stemness, neuronal hyperexcitability, and immunosuppression markers. Both agents potentiated the cytotoxic and apoptotic effects of temozolomide, supporting glutamatergic inhibition as a strategy for improving chemosensitivity. Conclusion: By establishing a framework for decoding glioblastoma heterogeneity through epigenetic aging, we identified the glutamatergic pathway as a clinically actionable vulnerability. Our findings suggest that combining anti-glutamatergic therapies with temozolomide exerts synergistic antitumor effects while mitigating adverse chemotherapy-induced phenotypes, such as increased stemness, neuronal hyperexcitability, and immunosuppression, thereby laying the groundwork for novel therapeutic strategies.

Immune checkpoint inhibitors targeting the PD-1/PD-L1 axis have shown great promise in treating bladder cancer and are now part of the standard treatment for advanced disease. However, many patients still fail to respond to treatment and at present many biomarkers are assessed but have yet shown only limited results. Therefore, with the advent of combination treatments and the increase of immune related adverse event, the search for reliable predictive biomarkers is paramount. Using a multiplexed enhanced FRET-FLIM based technique (QF-Pro) we quantified the interaction of PD-1/PD-L1, CTLA-4/CD80 and TIGIT/CD155 immune checkpoints in a pre-treatment TMA of 46 patients treated with atezolizumab. The association between higher PD-1/PD-L1 ICP interaction state and treatment efficacy was demonstrated in the male sample cohort, where it identified patients with better PFS. Conversely, patients exhibiting higher CTLA-4/CD80 engagement had a worse response to atezolizumab. Remarkably, the dual assessment of patients with high PD-1/PD-L1 and low CTLA-4/CD80 allowed to identify the best responders. These results indicate that the monitoring of patients immune profile in urothelial carcinoma might be critical in identifying patients who may benefit from combination therapy.

Liver metastases (LM) remain a major cause of death from different cancer types, in particular malignancies of the gastrointestinal tract. Liver metastases predict a poor response to immunotherapy due, among others, to the immunotolerant microenvironment (ME) of the liver and loss of local and systemic cytotoxic T cells. Thus, strategies that can reprogram the immune ME of the liver and restore cytotoxic T cell reactivity are being sought. We previously reported that estrogen signaling blockade impedes the growth of LM by reducing MDSC accumulation and monocyte/macrophage polarization. The aim of this study was to elucidate the underlying mechanism(s) and assess whether estrogen signaling in the myeloid lineage was driving the immunotolerant ME of LM. To this end, we generated mice with conditional myeloid cell-specific deletions of estrogen receptors (ER) or ER{beta} and analyzed in these mice the effect of ER loss on the liver immune ME and the outgrowth of LM. In mice with ER, but not with ER{beta} deletion, we observed a marked reduction in the growth of murine colon carcinoma MC-38 liver metastases as compared to their respective controls. Flow cytometry and immunohistochemistry revealed a decrease in macrophages that were polarized to the pro-tumorigenic M2-like phenotype and a concomitant increase in activated CD8+ T and NK cells relative to controls. Bulk RNAseq analysis performed on hepatic immune cells infiltrating the liver revealed changes in the expression of key cytokines/chemokines mediating immune cell recruitment, activation and polarization, including Ccl5 (upregulated) and Csf1 (downregulated). Taken together, the data suggest that ER signaling in myeloid-derived cells programs the immune landscape and contributes to an immunosuppressive and metastases-growth permissive ME in the liver.

Background Conversion surgery for stage IV gastric cancer is supported by a pooled overall survival hazard ratio of 0.36 (95% confidence interval 0.32-0.40) and, in the largest international cohort, median survival of 36.7 versus 12.5-13.8 months on chemotherapy. Survival is measured from diagnosis; the median diagnosis-to-gastrectomy interval is 124 days, which patients must survive to be counted surgical. Methods We simulated cohorts of 3,177 stage IV gastric cancer patients from published parameters: background median survival 14.5 months; median diagnosis-to-surgery interval 124 days (category-specific 92-174 days). Surgery had no effect (true hazard ratio 1.00 by construction). Data were analysed as the literature analyses them (exposure fixed at baseline, follow-up from diagnosis), and by time-varying Cox and landmark analysis. Confounding by indication was added in a second scenario. Results Under immortal time bias alone the naive analysis returned a hazard ratio of 0.794 (95% simulation interval 0.743-0.851), median survival 16.8 versus 12.8 months. Time-varying Cox recovered 1.000 and landmark analysis 1.000-1.004. Bias scaled with the interval: 0.849 at 92 days, 0.715 at 174 days. Adding confounding, the naive estimate fell to 0.601 (0.560-0.644) at strength 0.5 and 0.356 (0.323-0.385) at strength 1.5, overlapping the published estimate; median survival 21.9 versus 8.7 months. Correcting immortal time alone left residual bias (hazard ratio 0.439). Conclusions The reported survival advantage of conversion surgery is reproducible where the operation does nothing; published estimates cannot distinguish benefit from bias. Resolving this requires individual patient data analysed with methods that assign person-time correctly, or completion of JCOG2301.

Pancreatic ductal adenocarcinoma (PDAC) is an aggressive malignancy characterized by high metastatic dissemination, therapy resistance, and poor clinical outcome. Inhibitor of differentiation 1 or ID1, is frequently overexpressed in PDAC and is associated with tumour progression and adverse clinical outcome. However, the mechanisms governing its post-transcriptional regulation remain insufficiently characterized. Here, we identify tumour-suppressive miR-615-5p as a regulator of ID1 expression in PDAC. Integrative in-silico target prediction prioritized miR-615-5p based on seed complementarity and thermodynamic stability with the ID1 3' -UTR. Expression analysis of available PDAC clinical datasets revealed reduced miR-615-5p expression associated with increased ID1 expression. Direct association was validated using luciferase reporter assays, where miR-615-5p suppressed 3' -UTR reporter activity of ID1 in a sequence dependent manner, while mutation of the predicted binding site attenuated this effect. Further biotinylated-RIP and AGO2-RIP assays demonstrated the co-enrichment of ID1 transcripts and miR-615-5p with AGO2 associated RISC complexes, while AntimiR mediated inhibition of miR-615-5p perturbs association between miR/ID1 to AGO2, supporting interaction specificity. Functionally, modulation of miR-615-5p altered ID1 expression and impacted PDAC cell migration in vitro. Mechanistic analyses further indicated that the miR-615-5p/ID1 axis influences autophagic flux where miR-615-5p mediated inhibition of autophagy suppresses ID1 dependent cellular migration. Collectively, these findings define a previously uncharacterized miRNA-dependent regulation of ID1 expression and link this axis to autophagy-associated migratory responses in PDAC cells. The study expands the post-transcriptional regulatory landscape of ID1 and provides a possible mechanism where suppression of miR-615-5p leads to ID1 overexpression and subsequent poor clinical outcome in PDAC cells.

Myeloproliferative neoplasms (MPN) are chronic hematologic malignancies characterized by clonal myeloid expansion, inflammation, oxidative stress, and progressive bone marrow (BM) remodeling that may culminate in fibrosis and secondary acute leukemia. Here, we evaluated the therapeutic efficacy and the underlying mechanisms of melatonin (MT) and liposomal melatonin (nano-MT) in preclinical MPN models. MT selectively inhibited clonogenic growth of patient-derived peripheral blood mononuclear cells and induced pluripotent stem cell-derived CD34 hematopoietic stem and progenitor cells in comparison to healthy controls. This effect was associated with increased apoptosis, reduced reactive oxygen species (ROS), and decreased glucose uptake, independently of MT receptor signaling. Transcriptomic profiling of primary MPN CD34 cells revealed suppression of MYC targets, G2M checkpoint signaling, ROS, and glycolysis pathways. In co-culture models, MT reduced stromal -smooth muscle actin and phosphorylated SMAD2/3, indicating inhibition of TGF-{beta}-driven mesenchymal stromal cell-to-myofibroblast formation. In tamoxifen-inducible SclCreER;JAK2V617F mice, nano-MT achieved efficient spleen and BM targeting. Therapeutically, nano-MT reduced erythrocytosis, myeloid progenitor expansion, and BM IL-1{beta} levels. Longitudinal micro-computed tomography and histological analyses demonstrated normalization of BM architecture, reduced osteosclerotic remodeling and splenomegaly, decreased reticulin deposition and megakaryocyte numbers. In a dose-escalation study, nano-MT restored erythrocyte, hematocrit, and platelet counts and normalized megakaryocyte-erythroid progenitors. Combination treatment with ruxolitinib further reduced leukocytosis, neutrophilia, and monocytosis. Collectively, these findings demonstrate that (nano-)MT attenuates MPN and BM remodeling by targeting metabolic, inflammatory, and fibrotic pathways. This study provides the first evidence for a therapeutic benefit of nano-MT in MPN and establishes a rationale for further translational evaluation.

Lysophosphatidic acid (LPA) is widely recognized as an extracellular lipid mediator; however, the functional significance of intracellularly produced LPA remains poorly understood. Here, we investigated the regulatory mechanism and functional role of a LPA-producing lysophospholipase D GDE4, also known as GDPD1, in prostate cancer cells. GDE4 expression is induced under ER stress conditions in a PERK-dependent manner and requires the transcription factor ATF3. Disruption of GDE4 expression resulted in altered intracellular levels of LPA and LPA precursor lysophosphatidylethanolamine, accompanied by reduced cell proliferation. RNA sequencing and subsequent validation identified a set of genes downregulated in GDE4-depleted cells. Pharmacological inhibition experiments indicated that peroxisome proliferator-activated receptor and {gamma} (PPAR and PPAR{gamma}) signaling pathways contribute to the regulation of these GDE4-dependent genes. Collectively, our findings suggest that GDE4-dependent lipid remodeling is associated with PPAR/{gamma}-mediated transcriptional regulation under ER stress conditions. These results provide a potential framework for understanding the link between intracellular lipid metabolism and stress-responsive gene regulation.

Chemoradiotherapy (CRT) is the standard-of-care therapy for many solid malignancies, yet predictive biomarkers of treatment response remain limited. We identified a germline single nucleotide polymorphism (SNP) in an intrinsically disordered region of the lysine demethylase KDM3C/JMJD1C (p.S464T) that is associated with CRT outcomes in locally advanced rectal cancers (LARC) and head and neck squamous cell carcinoma (LA-HNSCC). In silico modeling with AlphaFold predicted S464T substitution influenced interaction between phosphorylated KDM3C and RNF8 FHA domain. In cellular models, conversion of S464 to T464 increased sensitivity to DNA-damaging agents. S464T substitution impaired damage-induced MDC1-RAP80 signaling and downstream RAP80-BRCA1 colocalization. SNP carrying cells impaired DNA repair causing genotoxic stress that is associated with increased cGAS-cGAMP innate immune signaling and increased apoptosis. Population analyses with the SNP highlighted an increase incidence of UV-induced skin and other cancers, linking inherited variation in the chromatin regulatory gene KDM3C to genome instability, cancer risk, and therapeutic vulnerability.

ZNF217 is an established oncogenic transcription factor that promotes cancer progression and therapeutic resistance; however, the mechanisms regulating ZNF217 protein abundance remain poorly understood. Here, we identify ubiquitin-specific peptidase 15 (USP15) as a critical regulator of ZNF217 stability and define a reciprocal USP15-ZNF217 signaling loop that sustains malignant phenotypes in ovarian cancer. Stable overexpression of ZNF217 in OVCA420 ovarian cancer cells enhanced proliferation, epithelial-mesenchymal transition, migration, invasion, and extracellular matrix adhesion. Notably, ZNF217 overexpression increased USP15 protein abundance without altering USP15 mRNA levels, whereas ZNF217 depletion reduced USP15 protein levels, suggesting post-transcriptional regulation. Conversely, USP15 depletion markedly reduced ZNF217 protein abundance while increasing ZNF217 mRNA levels, indicating that USP15 regulates ZNF217 predominantly at the post-transcriptional level. Proteasome inhibition restored ZNF217 protein levels following USP15 depletion, further demonstrating that USP15 promotes ZNF217 protein stability. Functionally, USP15 depletion in ZNF217-overexpressing ovarian cancer cells suppressed proliferation and multiple metastatic phenotypes, including migration, invasion, extracellular matrix adhesion, anoikis resistance, and multicellular aggregate formation. In vivo, USP15 depletion significantly reduced tumor progression and metastatic burden and prolonged survival in mice bearing ZNF217-driven ovarian tumors. Furthermore, USP15 depletion enhanced the sensitivity of ZNF217-overexpressing cells to carboplatin, paclitaxel, and doxorubicin. Collectively, these findings identify USP15 as an upstream regulator of ZNF217 protein stability and reveal a positive-feedback loop between USP15 and ZNF217 that reinforces oncogenic signaling. Targeting USP15 may therefore represent an indirect therapeutic strategy for suppressing ZNF217-driven ovarian cancer, particularly given the challenges associated with directly targeting oncogenic transcription factors.

Next-generation sequencing (NGS) for the detection of somatic variants has become the method of choice in a variety of molecular oncology fields and in the clinic. Its use ranges from sequencing entire tumor genomes and transcriptomes to targeted clinical diagnostic gene panels. The NYU Langone Genome PACT (Profiling of Actionable Cancer Targets, LG-PACT) assay is a qualitative in vitro diagnostic test that uses targeted next generation sequencing (NGS) of formalin-fixed paraffin-embedded (FFPE) tumor tissue matched with normal specimens from patients to detect gene alterations in a targeted panel covering 606 genes and the TERT promoter. Indications for testing are cancer (solid tumors and hematological malignancies) where a mutational profile from multiple genes would be informative for disease stratification, prognosis, or treatment options including targeted therapies and eligibility for clinical trials. The test is intended to provide information on somatic mutations including point mutations, small insertions/deletions (indels), and copy number aberrations for diagnostic and treatment decisions. LG-PACT is a United States Food and Drug Administration (FDA) cleared diagnostic test (510K: K202304). The clinical interpretation of sequencing data of molecular tumor markers from NGS encompasses automated variant calling tools with human interpretation. This final mostly manual review of data step is intensive, involving highly trained scientists, encompassing literature review, interpretation and clinical tier classification by pathologists, who then provide a complete molecular diagnostic report to the treating oncologists. We provide analysis of 1339 clinical genomic profiles from 31 different cancers and their subtypes, comprising of central nervous system (CNS) 792 (59%) cases (incl. meningioma, glioma and glioblastoma), with 267 (20%) cases predominantly of lung, pancreatic and colorectal and 280 of others (21%). Here, we present the technical challenges of validating an NGS oncological diagnostic targeted assay for clinical grade accuracy and sensitivity for patient care. We show how copy number alterations provide a more comprehensive description of the tumors genomic profile. We then outline the utility of targeted panel sequencing based on certified pathologist selection of reportable variants for our current patient cohort. Where analysis of variant detection has led to 49.4% (661/1339) of our clinical tumor samples containing mutations in known therapy targeted genes, 35.6% (477/1339) with mutation detected in other genes, and 15% (201/1339) cases being negative.

Background: Immunochemotherapy (ICT) is considered standard in regards to care for small-cell lung cancer (SCLC) in extensive stages, yet reliable biomarkers for treatment response remain elusive. While previous univariate analyses suggest specific peripheral lymphocyte subsets correlate with survival, the systemic immune response involves complex, multivariate interactions that require advanced analytical approaches. Methods: This paper analysed high-dimensional flow cytometry data from 32 patients with stage IV SCLC treated with carboplatin, etoposide, and atezolizumab. Peripheral blood was analysed at baseline (V0) and longitudinally during treatment. To identify potential early predictive biomarkers and mitigate sample attrition in later cycles, we focused on baseline and measurements after two cycles of ICT (V1). We employed a rigorous machine learning framework utilising nested cross-validation, bootstrapping, and permutation-based statistical testing to evaluate eleven different regression and survival models. Results: Under model-appropriate metrics, regressors did not generalise (R2 <0); conversely, censoring-aware Random Survival Forests (RSF) successfully extracted robust prognostic signatures. Baseline immune profiles (V0) achieved a concordance index (C-index) of 0.66 (p= 0.015), while dynamic changes from V0 to V1 ({triangleup}V) achieved a C-index of 0.65 (p= 0.022). Crucially, absolute values measured after two cycles of ICT (V1) yielded no significant signal (p= 0.445). Feature importance analysis confirmed the prognostic value of Th17 normalisation and identified Naive Regulatory T cells and Memory B cells as candidate components. Conclusion: Machine learning validation confirms a predictive signal in the peripheral immune profile of SCLC patients. Early dynamic shifts in the balance between regulatory and effector immune arms are associated with prognosis, contrasting with the lack of signal in absolute counts after two cycles of ICT. These findings establish a proof of concept for multivariate liquid biopsy immune profiling, warranting confirmation in larger cohorts and highlighting the necessity of integrating systemic and tumour-intrinsic data.

MYCN amplification drives replication stress in high-risk neuroblastoma, yet how MYCN-amplified tumour cells tolerate this stress to sustain proliferation remains poorly understood. Here we show that FBXL12, an SCF ubiquitin ligase substrate receptor that targets the Fanconi anaemia protein FANCD2 for degradation at replication forks, as well as the broader Fanconi anaemia and replication stress transcriptional program are elevated in high-risk and MYCN-amplified neuroblastoma. High FBXL12 expression independently predicts poor survival across neuroblastoma patient cohorts. FBXL12 loss stabilizes FANCD2 on chromatin, elevates ATR-dependent replication stress signalling and DNA damage during S phase, and impairs proliferation of MYCN-amplified neuroblastoma cells in vitro and in vivo. Mechanistically, MYCN directly engages the FBXL12-FANCD2 complex and antagonises FBXL12-mediated degradation of FANCD2 at replication forks, revealing that the oncogenic driver of replication stress also actively preserves the chromatin-bound FANCD2 pool required to tolerate it. Beyond S phase, FBXL12 loss disrupts FANCD2-dependent mitotic DNA synthesis and transmits unresolved replication intermediates into daughter cells. FBXL12-deficient cells consequently show transcriptional activation of MYC target gene, ATR, and mTOR signalling programs, and this pathway-concordant state confers differential sensitivity to ATR, and mTOR-targeting compounds, nominating candidate therapeutic strategies for this disease subset. Together, these findings define a MYCN-FBXL12-FANCD2 axis as a clinically relevant vulnerability in high-risk neuroblastoma.

Bladder carcinoma in situ (CIS) is a multifocal, non-muscle-invasive disease with a high risk of progression to muscle-invasive cancer. Current management strategies are often guided by genomic profiling of single tumor samples, which incompletely capture tumor heterogeneity and may contribute to treatment failure. In particular, the multifocal nature of CIS raises uncertainty regarding the uniformity of genomic, immunologic, and microenvironmental features across anatomically distinct sites within the same patient. To address this, we performed spatial transcriptomic profiling of CIS-containing tissue from four anatomically distinct sites within a single individual. Unsupervised clustering with marker-based annotation, integrated with metabolic inference, identified epithelial tumor populations alongside stromal, immune, and smooth muscle compartments. While key cellular states were conserved, their spatial organization and relative abundance varied by site. Metabolic analysis further revealed region-specific microenvironments shaped by local cellular architecture. These findings indicate that both cellular composition and metabolic activity are spatially structured. Collectively, these results demonstrate that CIS exhibits significant intra-patient heterogeneity not captured by single-site profiling. These findings require validation in larger cohorts but support multi-region sampling could help improve risk stratification, biomarker development, and prediction of response to intravesical therapies, with potential implications for more personalized treatment strategies.

Pancreatic ductal adenocarcinoma (PDAC) harbors a distinct intratumoral microbiome. Yet rigorous characterization of its composition is hampered by pervasive environmental and procedural contamination. Sources of contamination have not been thoroughly explored, and the methods of decontamination have not been sufficiently evaluated in a benchmarking manner. We systematically collected >300 negative control (NCT) samples comprising paraffin from formalin-fixed paraffin-embedded (FFPE) samples, lysis buffer and sterile water over a period of four years processed by different laboratory persons (LP). All samples were sequenced using full-length 16S rRNA gene sequencing with Oxford-Nanopore Technologies. We benchmarked four decontamination methods (restrictive filtering, decontam, SCRuB, and the Nejman et al.-derived (Nj) pipeline) against fresh-frozen tumor samples (FF) from LSL-KrasG12D/+;LSL-Trp53R172H/+;Pdx-1-Cre (KPC) mice, using the abovementioned contamination assessment to calculate a composite score for the assessment. Further, we validated those methods via technical replicates. Microbial profiles of NCT samples were significantly determined by control type, LP, year and season reflecting complex batch effects. The 15 most abundant contaminants spanned well-characterized environmental taxa and human commensals from the oral cavity. The LP processing samples left a significant microbial trace highly contributing to the batch effect. Decontamination benchmarking demonstrated that the Nj method consistently outperformed alternatives in both composite score and inter-replicate concordance. Application of Nj to fresh frozen PDAC samples substantially reduced contaminant burden while preserving putative tumor-associated signals in FF but not FFPE samples. Our results support the adoption of the Nj decontamination approach for future intratumoral microbiome studies in fresh frozen tumor samples.